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ABSTRACT 

We explore the relative importance of the stellar mass density as compared to the inner dark halo, 
using the observed gas kinematics throughout the disk of the spiral galaxy NGC 4254 (Messier 99). We 
perform hydrodynamical simulations of the gas flow for a sequence of gravitational potentials in which we 
vary the stellar disk contribution to the total potential. This stellar portion of the potential was derived 
empirically from color corrected K-band photometry reflecting the spiral arms in the stellar mass, while 
the halo was modelled as an isothermal sphere. The simulated gas density and the gas velocity field are 
then compared to the observed stellar spiral arm morphology and to the Ha gas kinematics. We find 
that this method is a powerful tool to determine the corotation radius of the spiral pattern and that it 
can be used to place an upper limit on the mass of the stellar disk. For the case of the galaxy NGC 
4254 we find Rcr = 7.5 ± 1.1 kpc, or Rcr = 2.1R cxp (K'). We also demonstrate that for a maximal 
disk the prominent spiral arms of the stellar component over-predict the non-circular gas motions unless 
an axisymmetric dark halo component contributes significantly (> 1/3) to the total potential inside 2.2 
K-band exponential disk scale lengths. 

Subject headings: galaxies: kinematics and dynamics — galaxies: structure — galaxies: halos — 
galaxies: spiral — galaxies: individual (NGC 4254) 



1. INTRODUCTION 

In almost all galaxy formation scenarios non-baryonic 
dark matter plays an important role. Today's numerical 
simulations of cosmological structure evolution reproduce 
fairly well the observed distribution of galaxy properties 
in the universe (e.g. Kauffmann et al. 1999) and attempts 
to model the formation of single galaxies have been made 
as well (Steinmetz and Mullcr 1995). In these simulations 
the baryonic matter cools and settles in the center of dark 
halos where it forms stars. The distribution of stars and 
gas in a galaxy depends strongly on the local star forma- 
tion and merging history. At the same time that the stars 
are forming the halos evolve and merge as well. 

The final relative distribution of luminous and dark mat- 
ter in the centers of the resulting galaxies is under debate 
because the mass distribution of the dark matter compo- 
nent is difficult to assess directly. Measuring luminous 
and dark matter mass profiles separately requires inno- 
vative strategies because the halo is poorly constrained 
and equally good fits to measured rotation curves can 
be achieved for a wide range of visible mass components 
(Broeils and Courteau 1997). In order to define a unique 
solution to this so called 'disk- halo degeneracy', the 'max- 
imal disk' solution was introduced. It assumes the highest 
possible mass-to-light ratio (M/L) for the stellar disk (van 
Albada et al. 1985, 1986). A practical definition is given 
by Sackett (1997) who attributes the term 'maximal' to 
a stellar disk if it accounts for 85% ± 10% of the total 
rotational support of the galaxy at R = 2.2 R oxp . This 
approach has proven to be very successful in matching ob- 
served H I and Ha rotation curves (van Albada et al. 1985; 
Kent 1986; Broeils and Courteau 1997; Salucci and Persic 
1999) and also satisfies some dynamical constraints, such 
as the criteria of forming m = 2 spirals (Athanassoula et 
al. 1987) as well as observational constraints on the struc- 
ture of the Milky Way (Sackett 1997). However, modern 



numerical N-body simulations find significant central dark 
matter density cusps (Fukushige and Makino 1997; Moore 
et al. 1999). Even if the prediction of these strong den- 
sity cusps may not be entirely correct, the simulations find 
that the dark matter is of comparable importance in the in- 
ner parts of galaxies (Blumcnthal et al. 1986; Moore 1994; 
Navarro et al. 1996, 1997) and it thus has a considerable 
influence on the kinematics. In this case a stellar disk of 
a galaxy would turn out to be 'sub-maximal'. 

It is important to determine the relative proportion of 
dark and luminous matter in galaxies for a better under- 
standing of the importance of the baryonic mass in the 
universe. This proportion also bears information on the 
dynamics and structure of the dark matter itself. Spi- 
ral galaxies are well suited to study dark matter distribu- 
tions because their distinctly ordered kinematics provide 
an excellent tracer of the gravitational potential in the disk 
plane. Since bars in galaxies are very prominent features 
with distinct dynamic characteristics, they are especially 
well suited to evaluate the amount of luminous matter. So- 
phisticated studies of barred galaxies indicate that their 
stellar disks alone dominate the kinematics of the inner 
regions - the stellar contribution is 'maximal' (Debattista 
and Sellwood 1998, 2000; Weiner et al. 2001b). However, 
studies of our own Milky Way, maybe also a barred spi- 
ral, still do not give a clear answer as to whether the disk 
is maximal (Sackett 1997; Englmaier and Gerhard 1999) 
or not (Kuijken 1995; Dehnen and Binney 1998). Bot- 
tema's analysis of the stellar velocity dispersion in various 
galactic disks led to the conclusion that disks cannot com- 
prise most of the mass inside the radial range of a few 
exponential scale lengths (Bottema 1997). Aside from the 
dynamical analysis of single systems, other attempts to 
tackle this problem have been undertaken. Mailer et al. 
(2000) used the geometry of gravitational lens systems to 
probe the potential of a lensing galaxy. They concluded 
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Fig. 1. — K'-band Image of NGC 4254 with a total exposure time 
of 20 minutes at the Calar Alto 3.5m telescope. Bright foreground 
stars are masked out. The overlay shows the slit orientations of the 
spectrograph. We took 8 longslit spectra (angles labelled in bold 
font) crossing the galaxy's center to measure the 2D velocity field. 



that a maximum disk solution is highly unlikely. Courtcau 
and Rix (1999) applied statistical methods to learn about 
the mass distribution in galaxies. In their analysis they 
found no dependence of the maximum rotation velocity 
on a galaxy's disk size. This is considered to be a strong 
argument to rule out a maximum disk solution. The con- 
flicting findings of different studies leave the question of 
the relative proportion of dark and luminous matter in 
galaxies still open. 

In this paper we want to exploit the fact that the stellar 
mass in disk galaxies is often organized in spiral arms, thus 
in kincmatically cold non-axisymmetric structures. In the 
canonical CDM cosmology the dark matter is collision-less 
and dominated by random motions. Although the intro- 
duction of weakly self-interacting dark matter was pro- 
posed to avoid current shortcomings of the CDM model 
(Spcrgel and Steinhardt 2000) it seems to raise other, com- 
parably severe problems (Miralda-Escude 2000; Ostrikcr 
2000). Hence it seems reasonable to assume that CDM 
is not substantially self-interacting, but dynamically hot 
and thus not susceptible to spiral and non-axisymmetric 
structure. 

In light of this, the key to measuring the baryonic to 
dark matter fraction is to make use of the non-axisymmetric 
structure that we can observe in the stellar light distribu- 
tion. Using deviations from axisymmetry of stellar disks, 
several efforts have already been made to constrain the 
dark matter content of galaxies (eg., Visser 1980; Quillen 
1999, and references therein). Some of the most significant 
conclusions came from studies of massive bars, which are 
the strongest non-axisymmetric structures in disk galax- 
ies. Spiral arms comprise a less prominent, but still signif- 
icant mass concentration. Already very early theoretical 



calculations of gas shocking in the gravitational potential 
of a spiral galaxy (e.g. Roberts 1969) told us that we 
could expect "velocity wiggles" with an amplitude of 10 
to 30kms~ while crossing massive spiral arms. For ion- 
ized gas, measurements of the velocity to this precision can 
be achieved with common longslit spectrographs. The im- 
print of the spiral pattern in the velocity field of real galax- 
ies is indeed not very strong, as apparent in the 2D velocity 
fields of M100 (Canzian and Allen 1997), of low surface 
brightness galaxies (Quillen and Pickering 1997) or of a 
sample of spiral galaxies Sakamoto et al. (1999). There are 
only a few spiral galaxies without bars that show stronger 
wiggles in the velocity field that are associated with the 
arms, eg. M81 (Visser 1980; Adler and Westpfahl 1996) 
and M51 (Aalto et al. 1999). In order to still achieve the 
goal of measuring mass-to-light ratios we need to compare 
the weak features that we expect in the measured velocity 
field to detailed kinematic models. Using new high reso- 
lution K-band photometry to map the stellar component 
and employing a modern hydro-code to simulate galactic 
gas flows we are confident that our models show enough 
details and we can eventually measure mass-to-light ratios. 
If the arms are a negligible mass concentration relative to 
the dark matter distribution in the galaxy, these wiggles 
should appear only very barely in the velocity field. The 
main aim of this project is to find out what fraction of the 
rotation speed comes from a mass component with spiral 
arms. In order to do this we compare the strength of the 
wiggles in a galaxy's observed velocity field to a model of 
the gas velocity field arising in a potential whose disk-halo 
fraction is known. As input for our gas dynamical simu- 
lations, we derive the stellar potential of the galaxy from 
color-corrected K-band photometry, while the dark mat- 
ter component is modelled as an isothermal sphere with a 
core. Simulations are performed for a variety of potential 
combinations and values for the pattern speed of the spi- 
ral structure. The results from these simulations arc then 
compared to the observed kinematics. 

The paper is organized as follows. In the next Sections 
we present the galaxy sample and the observations (§2) 
and the general modelling techniques (§3). In Section 4 
we present results for the galaxy NGC 4254 (Messier 99) 
and discuss them in §5. We summarize in §6. 

2. THE SAMPLE DATA AND NGC 4254 

We aim to perform this study on a sample of about 
half a dozen galaxies to account for systematic errors in 
the analysis of single systems. We could not identify a 
data set in the literature that fulfilled satisfactorily the 
requirements for this project. Hence we decided to collect 
data on a sample of galaxies ourselves. The data acquisi- 
tion procedure was carried out in the following way. First 
we took NIR photometry of possible sample candidates to 
check whether their spiral structure shows up reasonably 
well in the K-band. From that sample we selected the 
most promising candidates and obtained kinematic data 
for them on a second observing run. 

The requirement for the kinematic measurements was 
to trace the gas velocity perturbations of the spiral arms, 
ideally across the whole disk. The two classical methods 
for obtaining 2D gas velocity fields are H I or CO radio ob- 
servations and Fabry-Perot interferometry. Alternatively, 
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longslit spectra taken at different position angles can be 
used to map the disk. 

Single dish H I or CO observations are not suited for 
our project because they suffer from relatively bad angular 
resolution and poor sensitivity to faint emission between 
spiral arms and the outer part of the disk. Sakamoto et al. 
(1999) recently published CO observations of a sample of 
spiral galaxies. In the velocity maps the signature of the 
spiral arms is in the majority of cases not or only barely 
visible. This most likely results from beam smearing of 
the anyway weak velocity perturbations. 2D Fabry-Perot 
velocity fields provide the required angular resolution but 
usually give only a very patchy representation of the disk: 
Mainly the H II regions show up in the map (e.g., Canzian 
et al. 1993; Regan ct al. 1996; Weiner et al. 2001a). Be- 
cause the coverage we can achieve by taking eight longslit 
spectra across a galaxy's disk is reasonably high, in com- 
bination with its good angular resolution and sensitivity 
to faint Ha-emission we chose to collect our kinematic in- 
formation by taking longslit spectra. 

This paper deals only with the results of the first galaxy 
from the sample: NGC 4254. 

2.1. Observations and Data Reduction 

Since we require photometric as well as kinematic data 
for this study, there are some constraints which apply for 
the sample selection. To reasonably resolve the structures 
of the stellar disk we prefer galaxies with low inclinations 
with respect to the line-of-sight (LOS). However, the LOS 
component of the circular motion increases with inclina- 
tion i. Due to the fact that the projection of the galaxy 
scales with cos(i) and the LOS fraction of the velocity 
scales with sin(i), galaxies in the inclination range between 
30° and 60° are best suited for yielding all the photometric 
and spectroscopic information. We also chose the sample 
to consist mostly of non-barred high luminosity galaxies 
with strong spiral arms. All data were acquired at the 
Calar Alto observatory's 3.5 m telescope. 

For all sample galaxies we decided to take near-infrared 
(NIR) photometry to study the luminous mass distribution 
and gas emission line (Ha) spectroscopy to acquire kine- 
matic information for the systems. The NIR images were 
taken during two observing runs in May 1999 and March 
2000 with the Omega Prime camera at the 3.5 m's prime 
focus (Bizenberger et al. 1998). It provides a field of view 
of 6.'76 x 6f76 with a resolution of 0'.'3961 per pixel. We 
used the K' filter which has a central wavelength of 2.12 
|im. The integration time was 20 minutes on target which 
allows us to trace disks out to about 4 scale lengths. Our 
observation sequence included equal amounts of time on 
the targets and sky frames used for background subtrac- 
tion. The data reduction was performed using standard 
procedures of the ESO-Midas data reduction package. In 
total, we collected NIR photometry data for 20 nearby spi- 
ral NGC-galaxies. From their NIR appearance we selected 
half of the galaxies as a subsample for which kinematic 
data should be taken. 

We obtained the gas kinematics from longslit spectroscopy 
measurements of the Ha emission. With the setup we 
used, the TWIN spectrograph achieved a spectral reso- 
lution of 0.54 A per detector pixel, which translates to 
24.8 km s~ LOS- velocity resolution per pixel, allowing us 



to determine LOS- velocities with ~ 7kms~ precision. 
We sampled the velocity field of the galaxy along 8 slit 
position angles 1 , all crossing the center of the galaxy (cf. 
Fig. 1). The slit of the TWIN spectrograph measures 
4'x l'/5 on the sky. The spectra span a region of 6130 
- 7170 A centered on Ha. The spectra were obtained 
during three observing runs in June 1999, May and De- 
cember 2000. The integration time for each single spec- 
trum varied between 600 s and 1800 s due to weather and 
scheduling constraints. For the data reduction we used 
the IRAF package and applied standard longslit reduction 
procedures. We determined the LOS velocity component 
of the ionized gas as a function of radius from the galac- 
tic center from doppler shifts of the Ha-line. The center 
of the emission line was determined by fitting separately 
a single Gaussian profile to it and the brighter N II-line 
at 6584 A. The weighted comparison of the two fits pro- 
vided us with the uncertainty of the line center position. 
Finally, each of the eight spectra were folded at the center 
of the galaxy to get the two rotation curves. To do this 
coherently, we determined an average wavelength distance 
D\ between a prominent sky line (at \ s ky) and Ha at our 
best guess for the galactic center from all eight spectra. 
Then, separately for each slit, the Ha rest wavelength in 
the center of the galaxy - where v c = - was assigned 
to be at X s ky + D\ . There the two rotation curves should 
be separated. This provided us finally with 16 rotation 
curves at different position angles per galaxy, each reach- 
ing out to about 2 arc minutes. All the 16 rotation curves 
are displayed later in Section 4, Fig. 6 together with the 
results from one simulation. 

2.2. Structural Properties of NGC 4254 

We selected NGC 4254 (Messier 99) as the first galaxy 
from our sample to be analyzed, because it shows a clear 
spiral structure with high arm- inter-arm contrast. NGC 
4254 is a bright Sc I galaxy located in the Virgo galaxy 
cluster with a recession velocity of 2407kms~ 1 adopted 
from NED 2 . We assume a distance of 20 Mpc towards 
NGC 4254, taken from the literature (Sandage and Tam- 
mann 1976; Pierce and Tully 1988; Federspiel et al. 1998). 
It has a total blue magnitude of Bt = 10.44 and a diam- 
eter of 5.4 x 4.7 arcmin on the sky. At 20 Mpc one arc 
second is 97 parsecs in the galaxy which translates to 38.4 
pc per detector pixel. Our Ha rotation curve for NGC 
4254 (the Ha rotation curve is plotted later in Fig. 4 to- 
gether with modelled rotation curves) rises steeply out to 
~ 35" (3.4 kpc) and then flattens at a rotation velocity 
of ~ 155kms _1 . This agrees well with earlier estimates 
(Phookun et al. 1993). From the kinematics we know that 
the south-western part of the galaxy is approaching. If we 
assume trailing spiral arms then the galaxy rotates clock- 
wise when viewed from our perspective. 

In the K-band NGC 4254 shows very prominent two- 
arm spiral features at most radii. The northern arm bi- 
furcates at R w 4.5 kpc, causing a three arm pattern in 
the outer disk. Furthermore, the galaxy exhibits consider- 
All position angles quoted in the text are in degrees eastward from 
north 

2 The NASA/IPAC Extragalactic Database (NED) is operated by 
the Jet Propulsion Laboratory, California Institute of Technology, 
under contract with the National Aeronautics and Space Adminis- 
tration. 



Probing for Dark Matter in Galaxy Disks 



c 

D 

15 



I 
O 

o 

CD 

cn 
D 

£ 





K- 


-b 


an 


d 


Brightness profi 


le 


0000 








1 




1 ; 


1000 


r%. 








original profile 


-- 


100 


— 










- 


10 


r 








corrected profile ^w 


^> 


1 


a 










1 . 



Azimuthal profile at radius R 



20 40 60 80 100 120 140 
radius in [arcsec] 




Fig. 2. — Different results from the point-by-point correction of the mass-to-light ratio in the K' image, (a) The original azimuthally 
averaged brightness profile together with the corrected profile. The outward blueing of the galaxy leads to a stellar disk scale length which 
is smaller than in the uncorrected K'-band light, (b) The azimuthal brightness profile at a radius of 39'/2; azimuth is to the north. The 
profile shows a clear two arm structure that is preserved after the correction. Note, however, that the shape of the profile changes. 



able lopsidedness. NGC 4254 shows a strong arm-interarm 
brightness contrast, noted by Schweizer (1976) to be even 
stronger than the one for NGC 5194 (Messier 51). Even in 
the K-band the brightness contrast is rather high, approx- 
imately a factor of 2 over a wide radial range. Gonzalez 
and Graham (1996) argued that, combined with a usual 
density wave, some external mechanism is needed to in- 
voke such high contrasts. However, NGC 4254 is well sep- 
arated from any other galaxy in the cluster and there are 
no obvious signs for recent interaction. Phookun et al. 
(1993) reported in a recent paper the detection of high 
velocity H I clouds outside the disk's H I emission. The 
authors argue that in-falling H I gas may be responsible 
for the unusual 'one armed' outer structure of the spiral. 
In that case the spiral structure of NGC 4254 may have 
been recently reorganized, enhancing the southern arm or 
the arm-interarm mass density distribution. Recent inter- 
actions could in principle corrupt the project's assumption 
of a steady state spiral pattern. But as we will find from 
our full analysis, the steady state assumption is seemingly 
not far off on time scales of a few dynamical periods. 

NGC 4254 harbors a small bar-like structure at its cen- 
ter with a major axis position angle of w 40°. From both 
ends of the bar two major arms emerge with a third arm 
splitting off the northern arm. By analyzing a g — K color 
map of NGC 4254 one learns that this third arm is signifi- 
cantly bluer than the other regions of the galaxy and thus 
consists of a younger stellar population. 

In the K-band the disk of NGC 4254 is well approxi- 
mated by an exponential with a scale length of R exp ~ 
36", corresponding to ~ 3.5 kpc if a brightness average of 
the arm and interarm regions is considered and bright H II 
regions are removed from the image. Except for the very 
center, the whole surface brightness profile is well fitted 



by a double exponential model with an inner 'bulge' scale 
length of w 0.6 kpc (cf. Fig. 2a) 

3. MODELLING 

In the forthcoming Sections we describe the modelling 
needed to connect our photometric and kinematic obser- 
vations. In this pilot study, we restrict our analysis and 
description to the spiral galaxy NGC 4254 (Messier 99). 
The data modelling involves four discrete steps. First we 
derive a stellar gravitational potential from the K-band 
images. Second, for each stellar mass-to-light ratio, we 
find a dark matter halo profile to match the total rotation 
curve. Third, we perform hydrodynamical simulations of 
the gas within the combined stellar and dark matter po- 
tential. Finally, we compare the predicted gas velocity 
field to our Ha-observations. 

3.1. The Stellar Potential 

3.1.1. K-band Light as Tracer of the Stellar Mass 

We would like to derive the stellar potential directly 
from the deprojected K-band stellar luminosity density of 
the galaxy. To this end, it is important to understand how 
well the NIR surface brightness traces the stellar mass den- 
sity. There are two major concerns which complicate the 
direct translation: both dust extinction and spatial pop- 
ulation differences could introduce K-band mass-to-light 
ratio variations. Observing at ~ 2 |j.m reduces the effect 
of the dust extinction significantly. Since we look at most 
galaxies from a nearly face-on perspective, we expect the 
optical depths to be ~ 0.5 in the K-band, leading to lo- 
cal flux attenuations of < 10% (Rix and Rieke 1993; Rix 
and Zaritsky 1995). Significant global M/L- variations may 
arise from red supergiants; they emit most of their light in 
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Fig. 3. — Different stages of processing of the K-band image prior to the calculation of the potential, (a) The original frame, (b) Image of 
the galaxy after the correction for constant mass-to-light ratio was applied, (c) Final image, after removal of compact star forming regions 
by means of azimuthal Fourier decomposition. The images are not yet deprojected to face-on. 



the near infrared and are fairly numerous and thus could 
considerably affect the total light distribution in K-band 
images. Supergiants evolve very rapidly and therefore are 
mainly found within the spiral arms where they were born. 
This would tend to increase the arm-interarm contrast, 
leading to stronger inferred spiral perturbations in the stel- 
lar potential. 

We used the approach of Bell and de Jong (2001) to 
correct for local mass-to-light ratio differences, which may 
still be present in the K-band images. From spiral galaxy 
evolution models, assuming a universal IMF, Bell and de 
Jong (2001) find that stellar mass-to-light ratios correlate 
tightly with galaxy colors. They provide color dependent 
correction factors, which we can use to adjust our K-band 
image to constant mass-to-light ratio. For NGC 4254 we 
used a g-filter image taken from Z. Frei's galaxy catalog 
(Frei et al. 1996). The ^-filter (Thuan and Gunn 1976) has 
a passband that is very similar to the Johnson ^-filter, so 
we can use the color correction given for V — K while using 
g — K colors. The mass-to-light ratio is given by 

log 10 (M/L) = -1.087 + 0.314(V - K). (1) 

This relative correction was performed on a pixel by pixel 
basis by scaling each flux by the corresponding M/L (cf. 
Fig. 3b). NGC 4254 has V-K of 2.61 which results, ac- 
cording to (1), in a K-band maximal-disk stellar M/L, 
Y*, of 0.54. If we estimate a maximal-disk stellar M/L 
from the disk potential we use for the simulations, we find 
Y* w 0.64. These values are in good agreement, consider- 
ing the uncertainties that contribute to the estimation of 
Y*. In our approach even small changes in NGC 4254's 
distance or inclination result in errors of w 20%. The 
galaxy becomes significantly bluer from the center to the 
outer disk regions, with the bluest regions being the spiral 
arms themselves. As a result the mass density profile of 
the corrected disk gets steeper: the disk scale length (after 
deprojection) changes from 3.5 kpc to 2.7 kpc (cf. Fig. 2a). 
In general the correction for constant mass-to-light ratio 
does not particularly change the arm-interarm contrast of 
NGC 4254 and thus does not reduce the non-axisymmetric 
component of the stellar potential. In the force field of the 
galaxy, the relative strength of the tangential component 
makes up ~ 6% of the radial force component before and 



after the correction. Nevertheless, the correction modi- 
fies the non-axisymmetric force component on average by 
~ 30% (cf. Fig. 2b). This has an evident effect on the 
data and the models, which we think is favorable. The 
correction seems to work fine for most of the disk. How- 
ever, there are some thin dust lanes at the inner edges of 
the spiral arms in the green image of NGC 4254 that are 
optically thick. Because of the quite poor resolution of the 
green image, the dust lanes are not distinctly visible in the 
image. For these the correction leads to some overcorrec- 
tion because virtually all the green light is absorbed by the 
dust lane. But since these optically thick dust lanes are 
narrow and still at the location of the NIR spiral arms, 
we assume that their influence on the total potential is 
negligible. 

To reduce the influence of star forming H Il-regions and 
OB-associations, which are distinctly visible in the K-band 
data, we filtered the image by means of a Fourier decom- 
position. From the azimuthal decomposition we discarded 
the Fourier terms higher than N=8 and subtracted the 
residual from the K-band image of NGC 4254. This cor- 
rection does not depend significantly on the number of 
Fourier terms included in the fit. This procedure cleaned 
the image of the patchy small scale structure and left us 
with the smooth global spiral pattern (cf. Fig. 3c). 

3.1.2. Deprojection 

To deproject the K-band images to face on, we need to 
derive the disk inclination, i, and the major axis position 
angle, PA from the observed gas kinematics. We perform 
a global x 2 " n t °f an axisymmetric rotation curve model to 
the 16 observed rotation curves. The model is based on a 
combination of stellar and dark halo rotation curves which 
were scaled to an average rotation curve from the six slit 
positions closest to the major axis (cf. Fig. 4). Deriving i 
and PA is an iterative process. The final values get deter- 
mined only after the first full hydrodynamical simulations, 
when actual non-axisymmetric model rotation curves can 
be fitted to the observed kinematics. As it turns out, in 
the case of NGC 4254 the difference in the resulting incli- 
nation can be considerable. An axisymmetric model yields 
a disk inclination of i = 30? 8 while the simulations suggest 
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Fig. 4. — Different decompositions of the total gravitational potential into stellar and dark halo components, illustrated for the example 
of NGC 4254. The measured rotation curve (data points with error bars) — an average of six single rotation curves closest to the major axis 
- is well fit by all decompositions (solid line — ). Displayed are three cases for the contribution of the stellar disk f<j = 20%, 60% and 100% 

(dash-dotted line — ■ — ). The halo contribution (dashed line ) was adjusted to yield a good fit to the measurements in combination with 

the pre-selected stellar contribution. The fit parameters for the three cases are listed in Table 1. 



an inclination of i = 41?2. We use the latter as the incli- 
nation of the galaxy. However, spatial deprojection effects 
are not particularly large in the inclination range of 30° - 
40°. The major axis position angle is not sensitive to the 
details of the rotation curve; we adopt PA = 67?5. 

3.1.3. Calculation of the Gravitational Potential 

The smoothed image is then used to calculate the max- 
imal disk stellar potential by numerical summation over 
the whole mass distribution of the galaxy, which is defined 
by the surface mass density £(r) we get from the K-band 
images. 



$*(r, 



3& 



(2) 



with: 



\rt -Tj\ = VAx 2 + Ay 2 + e 2 



(3) 

The indices i and j denote different radius vectors, speci- 
fying here the pixels of the image array. Beyond the size 
of the detector array, the surface mass density is assumed 
to be zero. The stellar mass-to-light ratio Y* is assumed 
to be constant, after the corrections described in §3.1.1. 
The factor e in Equation (3) accounts for the finite thick- 
ness of galactic disks following (Gnedin et al. 1995). This 
softening factor was chosen to be compatible with a verti- 
cal exponential density scale height of h z = 400 pc, inde- 
pendent of radius. This assumption most probably holds 
for NGC 4254, which is a late type spiral and its central 
spheroid, which appears barred and extends only to about 
a disk scale length. Additionally, small bulges may not be 
spherical and thus might be well described by the constant 
scale height of 400 pc assumed for most of the simulations. 
Moreover, our derived potential does not depend signifi- 
cantly on the choice of h z . Varying h z by 10% leads to 
relative changes of the resulting forces by ~ 5.1 x 10~ 3 on 
average. Only for the most massive stellar disk models we 
performed an explicit disk-bulge decomposition and sub- 
stantially increased the softening factor for the bulge to 
760 pc. This is helps to achieve a better fit at the very 
inner part of the observed rotation curve. 



We also calculate a potential $ iax from an equivalent 
axisymmetric density distribution fitted to the measured 
K-band luminosity profile. The rotation curve is evaluated 
from the stellar potential by 



<(R) = R 



d$*, ax (r) 



dr 



(4) 



The maximal stellar mass-to-light ratio in equation (2) is 
determined by fitting the rotation curve emerging from 
the axisymmetric potential to the inner part of the ob- 
served rotation curve with the highest possible Y*. Later, 
when combining stellar and dark halo potentials, we take 
fractions f<j of this stellar 'maximum disk' fit, to explore 
various luminous-to-dark matter ratios. 

3.1.4. The Dark Halo Potential 

For the present analysis we use the radial density profile 
of a pseudo-isothermal sphere as the dark matter compo- 
nent in our model. Its radial density profile, characterized 
by a core of radius R c and a central density p c , is given 
by: 

p{r)=pc i+(-^r 

tie. 



(5) 



The corresponding rotation curve is 

R 



1 



(Kent 1986) and the potential is 

$halo(i?) 



. R 

arctan 
R V R, 



R vl lo (r) 



dr. 



(6) 



(7) 



The asymptotic circular velocity in this (infinite mass) 

halo Voo is related to the two other free parameters R c 
and p c by 

Voo = y/4nG Pc R*. (8) 

Thus, Wqo and R c uniquely specify all properties of the 
halo. For any fraction fj of the maximal stellar mass we 
can obtain the halo parameters from the best fit of the 
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combined stellar and halo rotation curve to the observed 
kinematics. 

In the final potential the two components are combined 
in the following way: 

$ tot (#|l d ) = f d **(H) + $halo(H|fd) (9) 

with fd ranging from to 1, and <&* denoting the stel- 
lar potential with maximal Y*. The contributions to the 
circular velocity add quadratically: 

v 2 tot (R, f d ) = f d vl(R) + v 2 alo (R, f d ). (10) 

For every f d we are able to find a dark matter profile, which 
complements the luminous matter rotation curve. % 2 -fits 
to the observed data are similar (cf. Table 1), necessitating 
this study. 

3.1.5. The Pattern Speed fl p and Corotation 

The spiral structure in galaxies is tightly related to the 
density wave theory that is used to describe the formation 
of these structures in spiral disks (for a review cf. , Athanas- 
soula 1984). If the spiral pattern evolves more slowly than 
the orbital time scale, resonances occur at certain radii in 
the galactic disk. We assume that the spiral pattern re- 
mains constant in a particular corotating inertial frame, 
specified by the pattern speed fi p . At the corotation ra- 
dius Rcr (fi P = V/Rcr where V is the circular rotation 
speed) everything corotates with the spiral pattern. In 
addition to this corotation resonance there are more radii, 
at which the radial oscillation frequency k and the circu- 
lar frequency £1 are commensurate. In linear density wave 
and modal theories, these so called inner and outer Lind- 
blad resonances should inclose the radial range for which 
a to = 2 spiral pattern may develop and persist over a 
longer period of time (Lin and Shu 1964). On the other 
hand, non-linear orbital models for open galaxies indicate 
that the symmetric, strong part of the stellar spiral might 
end at the inner 4:1 resonance (Patsis et al. 1991). Finding 
the radial locations of these resonances in real galaxies is 
a delicate endeavor and there are no widely agreed upon 
ways to determine them. 

For our analysis the pattern speed Q p or its correspond- 
ing corotation radius, is crucial since it modifies the forces 
acting in a non-rotating inertial frame 

F(R) = V<S> tot {R) + fl p R 2 . (11) 

Table 1 
Dark halo parameters. 



fd 


R c 


Woo 


X 2 /N 


(M D /M Dmax ) 


(kpc) 


(kms -1 ) 




0.2 


1.08 


155 


2.45 


0.4444 


2.00 


150 


3.42 


0.6 


3.03 


150 


3.30 


0.85 


5.68 


155 


2.00 


1.0 


7.30 


155 


1.84 



Note. - - Dark halo parameters used to gen- 
erate the potentials used for simulations. The 
\ 2 /N- values refer to an axisymmctric model. 



Therefore we need to perform a set of hydrodynamical sim- 
ulations at different O p , to find the best matching value 
and thereby to also determine the corotation radius. As it 
turns out, the simulations are well suited to find these val- 
ues reliably. However, to reproduce all the gravitationally 
induced features found in the observed kinematics, we sus- 
pect that the assumption of a constant pattern speed for 
the whole disk is only a crude approximation of the dynam- 
ical processes occurring in the disk of NGC 4254. The pat- 
tern rotation speed most probably depends slightly on the 
radius. A radius dependent pattern speed would eventu- 
ally drive evolution of the spiral pattern. Although spiral 
structure is most likely a transient feature in the galactic 
disk, it persists for several dynamical time scales (Don- 
ner and Thomasson 1994; Patsis and Kaufmann 1999) 3 . 
But, since for now we have no means to reliably determine 
the radial dependence of the pattern rotation speed, we 
choose to ignore the time evolution of the spiral pattern 
and assume a constant pattern speed for the simulations. 

3.2. Hydrodynamical Simulations 

The next step is to calculate the expected kinematics 
for a cold gas orbiting in <J> to t with a certain fi p . We 
model the gas flow with the BGK hydro-code (Prender- 
gast and Xu 1993; Slyz and Prendergast 1999). This is an 
Eulerian, grid-based hydrodynamics code which is derived 
from gas kinetic theory. The simulations are challenging 
because cold gas in galactic disks rotates with such high 
Mach numbers that if the flow is diverted from a circu- 
lar orbit by non-axisymmetric forces, the consequence is 
etched in the gas in patterns of shocks and rarefied regions. 
BGK is well suited for these simulations because as veri- 
fied by extensive tests on standard ID and 2D test cases of 
discontinuous non-equilibrium flow (for a review, see Xu 
1998), at shocks and contact discontinuities BGK behaves 
as well as the best high resolution codes and it gives better 
results at rarefaction waves because it naturally satisfies 
the entropy condition. 

To model the two-dimensional gas surface densities and 
velocity fields for NGC 4254 we carried out a set of sim- 
ulations on a 201 by 201 evenly spaced Cartesian grid. 
Our data for NGC 4254 extend out to a radius of about 
11.6 kpc, hence on a 201 by 201 grid this gives a resolution 
of about 1 16 pc per side of a grid cell, which is considerably 
higher than the effective force smoothing of 400 pc. 

Given that the gas surface mass density of the mod- 
elled galaxies is much lower than the density of the stellar 
disk and halo, we neglect the self-gravity of the gas, and 
compute its response to a fixed non-axisymmetric gravita- 
tional potential derived from the corrected NIR image of 
the galaxy with a 3-fold reduced resolution (cf. Sect. 3.1). 
We approximate the ISM as a monatomic gas. 

The gas temperature profile is taken to be uniform and 
by imposing an isothermal equation of state throughout 
the simulation we assume that the gas instantaneously 
cools to its initial temperature during each updating timestep. 
The initial density profile of the gas is exponential with a 
scale length equal to a third of the disk radius, namely 3.86 
kpc. We begin each simulation with the gas initially in in- 

3 In light of our project it is interesting to note that in both N-body 
simulations referenced here, sub-maximal disks were adopted in the 
initial conditions 
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viscid centrifugal equilibrium in the axisymmetric poten- 
tial given by <&* iax and <l>haio (sec Sect. 3.1). Following the 
initialization of the gas in centrifugal balance in an axisym- 
metric potential, we slowly turn on the non-axisymmetric 
potential at a linear rate computed by interpolating be- 
tween the final non-axisymmetric potential and the initial 
axisymmetric potential so that the potential is fully turned 
on by the time 40 sound crossing times of the code have 
passed. Here the sound crossing time is defined as the time 
it takes to traverse the length of the diagonal of a grid cell 
at sound speed. For an isothermal simulation with a sound 
speed of 10 km s~ the sound crossing time of a cell in our 
simulation is about 16 Myr, so that by 40 sound crossing 
times, the gas has evolved for 640 Myr which is about 1.4 
times the dynamical time of the galaxy measured at a ra- 
dius of 11.6 kpc. After the non-axisymmetric part of the 
gravitational potential has been fully turned on, we con- 
tinue to run the simulation for about another dynamical 
time. A thorough description of the technical details of 
the simulations is given in a companion paper (Slyz et al. 
(2001), in prep.). 

We ran a large set of simulations both to understand 
the power and limitation of our modelling in general and 
to match the observations. Simulations were performed 
for a total of five different fractions f d of the stellar disk: 
disk only, i.e. f d = 1, and f d = 0.85, 0.6, 0.4444 and 0.2, or 
accordingly f d is given in % from 20% to 100%. In all the 
cases the core radius and the asymptotic velocity of the 
pseudo-isothermal halo were adjusted to best match the 
averaged rotation curve, as summarized in Table 1. The 
variations in \ 2 /N between the low mass disks and the 
massive ones are mainly caused by the attempt to keep 
R c and Woo at physically reasonable values. The bump at 
20" to 40" (Fig. 4) is fitted better for the high mass disks, 
which reduces the overall \ 2 compared to the low mass 
disks. 

We have no secure prior knowledge of the spiral pattern 
speed il p . We determine it by assuming different values 
for ilp and then comparing a simulation to the data. For 
every simulation with a different stellar/dark halo combi- 
nation, we get slightly different values for the best match- 
ing f2 p or equivalently for the corotation radius -Rcr- We 
covered the complete range of reasonable Rcr, i-e. from 
about a disk scale length to well outside the disk. We even 
made simulations for the case of no spiral pattern rotation, 
Rcr -» oo. 

To test how the amplitude of the velocity perturbations 
depends on the responsiveness of the gas, we ran simula- 
tions at a variety of temperatures corresponding to sound 
speeds c s , of 10, 20, 30 and 40 km s _1 . In the following Sec- 
tions we discuss some of the results from the simulations. 

4. RESULTS FOR NGC 4254 

4.1. Simulated Gas Density 

Figure 5 shows eight views of the simulated gas density 
for different pattern speeds over-plotted as contours over 
the unsmoothed, deprojected, color corrected K-band im- 
age of NGC 4254. The simulated gas density follows an 
overall exponential profile with a scale length of w 4.2 kpc, 
comparable to the one of the disk itself. The contours in 
the Figure are chosen to highlight the density enhance- 
ments and locations of the gas shocks caused by the spiral 



arms. For almost all simulated cases the strong part of the 
galaxy's spiral structure lies well inside corotation, where 
the circular velocity is larger than the spiral pattern speed. 
The gas will thus enter the spiral arms from their inward 
facing side, producing the strongest shocks there. For a 
well matching simulation, we expect the shocks to be near 
the OB associations that trace the spiral arms. 

It is remarkable how well the overall morphology of 
NGC 4254 can be matched by the gas density simulations. 
Not only are the two major spiral arms clearly identifi- 
able in most simulations but the less prominent north- 
ern arm and the locations where the arms bifurcate are 
reproduced in some cases well. For fast pattern speeds 
(-Rcr = 5.4 — 7.58 kpc in Fig. 5) we find a strong shock in 
the northern part of the galaxy that cannot be correlated 
with any mass feature. We believe it develops because the 
potential close to the upper boundaries of the computa- 
tional grid is quite non-axisymmetric, and this leads to 
a spurious enhancement of a shock. The shock does not 
propagate into regions inside the corotation radius, and 
therefore we refrain from smoothing the potential. 

It is important to note that the simulations lead to a 
very stable gas density distribution that does not change 
much after the non-axisymmetric potential is fully turned 
on. When the contribution of the disk is increased in the 
combined potential, all spiral features get enhanced in the 
gas density but the galaxy morphology is essentially un- 
changed. For a more detailed discussion of these issues, 
please refer to the accompanying paper by Slyz et al., 
which is in preparation. With increasing pattern rotation 
(smaller corotation radii in Fig. 5) we find that the pre- 
dicted spiral arms become more and more tightly wound. 
For a comparison to the stellar spiral morphology we need 
to define some criteria to pick the right model. If the sit- 
uation in NGC 4254 is similar to NGC 4321, whose gas 
and dust distributions and their connection to star form- 
ing regions have been discussed in detail by Knapen and 
Bcckman (1996), then a good matching gas morphology is 
one where for radii smaller than the corotation radius, the 
shocks in the gas density lie on the inside of the stellar spi- 
rals. Shortly downstream from there, many star forming 
H II regions, triggered by gas compressions, should show 
up as patches in the arms. According to these criteria, the 
best matching morphology 

can be unambiguously identified to be produced by a 
simulation with Rcr ~ 7.6 kpc (Fig. 5 upper right panel), 
corresponding to a pattern speed of Q p « 20 km s~ kpc~ . 
Rcr. = 6.4 kpc and Rcr. — 8.3 kpc enclose the range of 
possible values. This corresponds to an uncertainty of ~ 
15% in the value of Rcr- 

Our results were compared to values of Rcr for NGC 4254 
from the literature, which were determined by different 
means. The results Rcr ~ 8.45 kpc (Elmegreen et al. 
1992) and R C r = 10.2 ± 0.8 kpc (Gonzalez and Gra- 
ham 1996), scaled to our distance assumptions, provide 
larger estimates than our findings. However, given the 
picture of non-linear orbital models, where one expects 
the strong part of the stellar spiral to end inside the inner 
4:1 resonance (Patsis et al. 1991) we find that our result 
of Rcr ~ 7.6 kpc is consistent with the galaxy's morphol- 
ogy. In short, these simulated gas densities provide us with 
an excellent tool to determine the pattern rotation speed 
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Fig. 5. — Simulation results of the gas density distribution overplotted in red contours on the deprojected K'-band image of NGC 4254 
(center). From the image an axisymmetric radial brightness profile has been subtracted to enhance the contrast of the spiral arms. The 
results are displayed for eight different assumptions for the pattern speed, respectively the corotation radius i?CR- The green circle marks 
corotation. The range goes from fast rotation i?CR = 5.4 kpc to slow pattern rotation i?CR = 15 kpc, lying mostly outside the image frame 
and finally for no pattern rotation. 



of the galaxy. Apart from requiring a constant global pat- 
tern rotation our approach is independent of an underlying 
spiral density wave model. The overall very good repre- 
sentation of the whole spiral structure by the simulated 
gas density makes us rather confident that the simulations 
render realistic processes affecting the gas. 

4.2. Simulated Gas Velocity Fields 



As another output of our simulations we get the two- 
dimensional velocity field of the gas. As is evident from 
a comparison to the gas density distribution, the velocity 
jumps are - as expected - at the locations where the den- 
sity map shows the shocks. They show up as areas of lower 
local circular velocity compared to the elsewhere rather 
smoothly varying gas velocity field. The velocity wiggles, 
as well as the density shocks themselves, have very tight 
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profiles and thus are even more narrow than the physical 
extent of the stellar arms. They have to be compared to 
the observed kinematics. 

4.2.1. The observed kinematics 

The rotation curves at the 16 slit positions from the ob- 
servations are shown in Fig. 6 as data points. It is appar- 
ent that the longslit spectra allow a good velocity coverage 
along the slits. Almost all rotation curves show contigu- 
ous data points out to a radius of > 1.'5. The spectra 
show a lot of wiggles on a small spatial scale. Jumps of 
< 30kms~ on a scale of ~ 5" are common. The very 
prominent jumps that we observed in slit positions 22?5 
and 225° clearly exceed the average wiggle sizes. Inside of 
about 0f3 the small bar influences the velocity field. The 
most prominent trace of the bar occurs at its minor axis 
at the slit positions of 135°/157?5 and 315°/337?5. 

The trace of kinematic features in the outer disk is not 
conspicuous in subsequent slits. The eastern part of the 
disk (slit positions 45° - 135°) displays a quite smooth ve- 
locity field, while the western part shows some large scale 
variations. Aside from the inter-arm region between the 
inner disk and the southern arm where a « 0.'5 wide de- 
pression is moving outward in subsequent slits (positions 
> 247? 5) no significant features are apparent in the outer 
disk. Unfortunately in this inter-arm region the S/N is not 
so good. Does that mean we do not see the trace of the 
arms in the velocity field, or is the single slit representa- 
tion of the 2D velocity field misleading and does not allow 
us to identify coherent features in adjacent slits? Clearly, 
the wiggles associated with spiral arms in NGC 4254 are 
not nearly as strong as in M81 (Visser 1980), thus their 
identification is harder. A CO map of NGC 4254's cen- 
ter (Sakamoto et al. 1999) shows also no coherent wiggles 
across spiral arms and we doubt that it would be much 
different on a Fabry-Perot image. Rather than being con- 
fused by the one dimensional nature of our rotation curve 
slices we believe that the spiral features in the velocity 
field arc intrinsically weak. 

4.2.2. Projection and Alignment 

To compare the simulated gas velocities to the observed 
data, we need to project the modelled velocity field ac- 
cording to the real orientation of the galaxy so that we 
can obtain the line-of-sight velocity along locations corre- 
sponding to the slit positions taken with the spectrograph. 
For this procedure we take the velocity components and 
transform them into velocity components parallel and or- 
thogonal to NGC 4254's major axis. Since the simulations 
yield truly 2D velocity fields, the component parallel to 
the major axis does not contribute to the line-of-sight ve- 
locity as it reflects only tangential motion. The orthogonal 
component is multiplied by the sine of the inclination of 
the disk, to account for the line-of-sight fraction of the 
velocities. 

We read the velocities along slits which correspond to 
the slit positions along which we took measurements with 
the spectrograph. The angular width of the grid cells of 
the simulation (l'/19) is comparable to the slit width of the 
spectrograph we used (l'-'5), which was also the average 
seeing conditions. 

The angles between the spectrograph slit orientations 
must be translated to angles in the plane of the galaxy 



to actually compare the same parts of the velocity fields. 
This translation makes use of the following relationship: 

tan(<#nt - PA) • cos(i) = tan(</? app - PA) (12) 

where y a pp is the apparent angle of the spectrograph slit 
across the galaxy on the sky and ip- m t is the correspond- 
ing intrinsic angle within the plane of the galaxy. PA is 
the position angle of the galaxy's major axis on the sky. 
All angles are in degrees measured eastward from north. 
Solving for ip- lnt gives: 

^tan(v? app -PA)" 



Vint = arctan 



PA. 



(13) 



\ cos(i) 

The foreshortening of the radial proportions due to pro- 
jection in the direction of the minor axis is determined 

by 



-Rproj = R ' \/cOS 2 (<£>int) 



+ (sim> int )-cos(i)) 2 . (14) 

Finally the detector pixel sizes of the TWIN, where the 
observed velocities come from, and Omega Prime camera, 
whose pixel scale is the reference grid for the simulations, 
must be adjusted to perfectly align with each other. 

4.2.3. Overall Fit Quality 

Figure 6 shows the 16 separate rotation curves with a 
corresponding simulated velocity field over-plotted. The 
simulation used here for comparison is the one for which 
the gas density distribution best-fit the K-band image (dis- 
played in Fig. 5). It has a corotation radius of Rqr = 
7.58 kpc - corresponding to a pattern speed of Q p = 
20kms~ 1 kpc _1 - and a stellar disk mass fraction fd of 
60%. A gas sound speed of c s = lOkms -1 was assumed 
here. 

An overall fit quality gets determined by a global % - 
comparison of the simulated velocity field to the actual 
observed velocity field along the 16 measured slit posi- 
tions. Since with the errors from the observations we find 
a x 2 /N > 1.0 we added an overall 9.5 km s _1 to provide an 
average fit with x 2 /-/V w 1.0. The x 2 -htting excludes the 
very central region hosting the small bar because the mod- 
elling is not intended to fit the central bar, which might 
have a different pattern speed. The total number of data 
points included in the x 2 -fitting is 1077. 

The general fit quality is governed by the effect that the 
projection of the simulated velocity field introduces. The 
good overall match indicates that we quite reliably found 
the right position angle and inclination for the galaxy. The 
simulated velocities align very well with the measured data 
points. In addition to the good overall match, the gen- 
eral rising or falling shape of the separate curves is also 
excellently reproduced by the simulations. The lopsided- 
ness of the galaxy is reflected in the shape of the rotation 
curves on the receding and approaching side of the disk. 
At the receding side (67?5) the rotation curve rises steeply 
and continues to rise out to 2', while the approaching side 
(247?5) rises less steeply but flattens out or even drops be- 
yond 0?7. These characteristics are closely reproduced by 
the models. 

A close inspection of the two profiles shows however that 
the overlap in the match of the simulated velocities with 
the measurements is not always satisfactory. The agree- 
ment of local features in the simulations and the measured 
data is sometimes very good and even occurs in subsequent 
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Fig. 6. — Simulation results of the gas velocity field in comparison to the observed rotation curves. Displayed are the measured rotation 
curves (data points), the axisymmetric model (thin line) and the rotation curves from the hydrodynamical simulations (thick line) for all 16 
slit position angles. The parameters for the simulation were: fd of 60% and i?CR = 7.58 kpc. There arc no error bars plotted for the data, 
but the errors can be estimated from the point-to-point scatter of the data. 



slit positions. However, there are also many locations, 
where the match is poor. This is particularly the case in 
the inner region of the galaxy, where the small bar dom- 
inates the kinematics. Both profiles show strong wiggles 
where the slit crosses the bar, especially close to the minor 
axis of NGC 4254's velocity field, which is also close to the 
minor axis of the bar itself. While the simulations show 
a rather symmetric imprint, the measurements exhibit a 
signature different from that, leading to a significant mis- 
match at several slit positions, e.g. 292?5 and 337?5. This 
might be caused by the bar, having a slightly different 



pattern speed. In the outer parts of the rotation curve we 
also find several wiggles in the observed data that have no 
correspondence to the wiggles in the simulations and vice 
versa. 

It is important to note that we do not expect to re- 
produce all the wiggles in the galaxy's rotation curves, 
since we are only modelling those created by the non- 
axisymmetric gravitational potential. The wiggles origi- 
nating for another reason - like expanding SN gas shells 
- are not considered by the simulations and thus do not 
show up in the resulting velocity field. 
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Fig. 7. — Comparison of four simulations for different fractions f,j of the stellar mass component. The simulation was done with a corotation 
radius of 7.58 kpc and a gas sound speed of 10 km s~ 1 . Displayed is one of the 16 slit positions. Clearly the predicted "wiggles" in the rotation 
curve grow much stronger for higher disk mass fractions. 



One very interesting thing to mention is the fact that 
even for moderate differences in the fit quality of the ax- 
isymmetric disk model (cf. Table 1), the fit quality of the 
non-axisymmetric simulated velocity field is rather inde- 
pendent of the initial axisymmctric model. And, more- 
over, the formally preferred axisymmetric maximum-disk 
decomposition (cf. Table 1) turns out to be the most un- 
favored model, once the simulations were performed. This 
implies that even if an axisymmetric model profile pro- 
vides a better fit to a measured rotation curve, it does 
not necessarily mean, that this combination provides the 
best fit when one considers the 2D non-axisymmetric gas 
evolution. 

4.2.4. Varying the Stellar to Dark Matter Ratio 

As already mentioned in Section 3.2, we performed sim- 
ulations for five stellar disk and dark halo combinations, 
as listed in Table 1. Since the non-axisymmetric perturba- 
tions are induced in the potential by the stellar contribu- 
tion, we expect the amplitude of the wiggles in the mod- 
elled rotation curves to depend significantly on the non- 
axisymmetric contribution of the stellar potential whereas 
we expect the radial distribution of the wiggles to be rather 
independent of the stellar mass fraction. As expected, in 
the simulations with the lightest disk, the wiggles look like 
modulations on the axisymmctric rotation curve. In the 
case of the maximum disk, the rotation curves are strongly 
non-axisymmetric (cf . Fig. 7) . To describe this character- 
istic more quantitatively, we learn from Fig. 8 that the 
amplitude of the deviations from axisymmetry increases 
linearly with the mass fraction of the stellar disk, which 
proves the general validity of the concept with which we 
try to approach this problem. 

The strongest velocity wiggles arising in the modelled 
velocity fields are the ones connected to the central bar- 
like feature. However, since here we are not interested in 
modelling the dynamics of the bar, we exclude this inner 
part from the analysis. 

If we do a formal x 2 comparison of the models with the 
observed data we find that for most of the sub-maximal 
disks the fit is considerably better than for heavy stellar 
disks. By formal, we mean that we use all data points for 
the x 2- fit, regardless of whether a certain part matches 
well or not with two exceptions: we exclude the very cen- 
tral part with the bar and we correct for the outer strong 
shock appearing in the fast rotating models (cf. Fig. 5). 
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Fig. 8. — Deviation of the simulation from axisymmetry. Dis- 
played is the average deviation of each radial simulation bin from 
the axisymmetric rotation curve. It rises linearly with the stellar 
disk mass contribution f<j. 



The result from this X 2 -fit is presented in Fig. 9. In 
all cases, fd = 100% gives the worst fit to the observed 
rotation curves. For the lower mass disks it is very hard 
to decide, whether a particular disk mass is preferred. For 
-Rcr = 6.4 kpc we find about the same % 2 for all non- 
maximal disk models. Since we can not reject data on 
a physical basis, we can state only a trend at this point. 
Thus our conclusion from this part of the analysis is that 
the disk is most likely less than 85% maximal. 

4.2.5. Varying the Gas Temperature 

This kind of analysis may be very sensitive to the tem- 
perature of the gas which is assumed for the modelling. 
Higher gas temperature corresponds to a higher cloud ve- 
locity dispersion c s and thus to a reduction in the response 
of the gas to any feature in the gravitational potential. 
Within the Milky Way Galaxy c s varies from « 6kms _1 
in the solar neighborhood to w 25kms _1 in the Galactic 
center (Englmaier and Gerhard 1999). For simulations of 
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a galactic disk with an isothermal equation of state, the 
most commonly used values for c s are 8 - 10kms~ , cor- 
responding to < 10 4 K in gas temperature (Englmaier and 
Gerhard 1999; Weiner et al. 2001b). In these simulations 
the authors make the statement that within the reasonable 
limits of c s = 5 — 30 km s~ the modelled gas flows across 
the primary shocks is not considerably affected. Only 
when modelling strong bars in galaxies the simulations 
might be dependent on the choice of the gas sound speed 
(Englmaier and Gerhard 1997). For the main set of sim- 
ulations presented in this paper, we chose c s = lOkms - , 
corresponding to a gas temperature of 7250 K. To probe 
the effect of the gas temperature, we performed simula- 
tions for four different gas sound speeds, c s , of 15, 20, 
30 and 40 kms~ . With a simulation performed at higher 
sound speed, we tested if our results at lower sound speeds 
were caused by too cold gas being too strongly affected by 
the non-axisymmetric part of the potential. Increasing 
the sound speed slightly broadened and shifted the shock 
fronts in our simulations, but not enough to change the 
conclusions of this paper. A more detailed technical dis- 
cussion on the simulated gas dynamics in NGC 4254 will 
be presented in a forthcoming paper by Slyz et al. 

5. DISCUSSION 

To learn about the relative role of stellar and dark mat- 
ter in the inner parts of spiral galaxies, we have extended 
the hydrodynamical simulations of galactic features from 
bars (Weiner et al. 2001b) to the regime of spiral arms. In 
this paper we have used color corrected K-band images to 
derive the stellar portion of the potential of NGC 4254 for 
the analysis. 

From the simulation of the gas density we find a very 
good morphological match of the gas shocks to the spiral 



arms in the galaxy; this aspect of the simulations deter- 
mines the pattern speed fl p of the galaxy to about 15% 
precision. 

A comparison of the observed and simulated kinematics 
has turned out to be challenging. Although the overall 
shapes of the different rotation curves were very well re- 
produced by the simulations, some small scale structure 
remains unmatched. The formal comparison of the gas 
velocity field to the observed Ha kinematics favors simu- 
lations with small disk mass fractions fd (cf. Fig. 9) and 
correspondingly small values for the stellar mass-to-light 
ratio. With the K-band mass-to-light ratio discussed in 
Section 3.1.1 our results yield an overall stellar M/L of 
f * S 0-5. We can estimate the relative mass fractions 
from their contributions to the total rotation velocity. At 
a radius of 2.2 uncorrected K'-band exponential disk scale 
lengths (2.2i? cxp w 79" or 7.7 kpc), the individual rota- 
tional support of the stellar and dark halo components for 
fd = 85% are v* = 125kms~ and fhaio = 86kms _ . If a 
total mass is estimated via 
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we find that Mh a io ~ 0.47 M* at R — 2.2i? oxp , or accord- 
ingly > 1/3 of the total mass inside R cxp is dark. Since our 
confidence limits are not very tight, we cannot use them 
to test other authors' findings in detail. Projects yielding 
results in favor of a sub-maximal stellar disk usually find a 
disk mass fraction less than our upper limit estimate. Bot- 
tema (1997) as well as Courteau and Rix (1999) conclude 
that the contribution of the stellar disk to the total rota- 
tion is v* ~ 0.6i>tot which translates to Mhaio ~ 0.6M to t- 
At the current state of the project we cannot exclude or 
confirm these findings. This issue is going to be discussed 
more thoroughly as soon as we have a few more examples 
analyzed. 

In the following Sections we will discuss the details that 
could cause deviations from a perfect match between the 
simulations and the measurements. 

5.1. Is the Concept Reasonable? 

If the self-gravity of the stellar mass in the disks of spi- 
ral galaxies plays an important role then undoubtedly the 
potential becomes non-axisymmetric. The trajectory of 
any kinematic tracer in the galaxy, such as the gas, should 
be affected by these potential modulations. But is the 
H II component of the gas the best choice for tracing the 
galaxy's potential? Analytic calculations of gas shocks in 
the gravitational potential of a spiral galaxy (e.g. Roberts 
1969) tell us that we should expect velocity wiggles with an 
amplitude of 10 to 30 kms~ while crossing massive spiral 
arms. However, kinematic feedback to the gas from re- 
gions of massive star formation, from expanding gas shells 
produced by supernova explosions and from other sources 
of turbulence, introduces small-scale random noise in the 
velocity fields. These fluctuations typically lie in the range 
of 10 to 15 kms -1 (Beauvais and Bothun 1999) and seem 
even higher in the case of NGC 4254. The kinematic small- 
scale noise could be increased if the dynamics of the bright- 
est H II regions is kincmatically decoupled from the global 
ionized gas distribution. To check this we over-plotted the 
Ha-intcnsity on the rotation curves to see if the star for- 
mation regions coincide with the strongest wiggles in the 
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rotation curves. There is, however, no discernible relation 
between the amplitude of a wiggle and the intensity of the 
H II region, indicating that the deviations are not confined 
to compact H II regions. 

As an alternative to observing the ionized phase of the 
hydrogen gas, one could consider using H I radio observa- 
tions. However, available H I data are limited by the larger 
size of the radio beam that smears out kinematic small- 
scale structures in the gas. In principle, stellar absorption 
spectra could also provide relevant kinematic information, 
but this approach has two disadvantages compared to an 
approach using gas kinematics. First, the acquisition of 
stellar absorption spectra with sufficient S/N would take 
a prohibitively large amount of telescope time. Second, 
stellar kinematics cannot be uniquely mapped to a given 
potential; there are different sets of orbits resulting in the 
same observed surface mass distribution and kinematics. 
Thus, despite its apparent shortcomings H II measure- 
ments seem to be the most promising method with which 
to approach the problem. 

Could the discrepancies between modelled and observed 
kinematics be related to taking the NIR K-band image of 
a galaxy as the basis to calculate its stellar potential? As 
already discussed in Section 3.1 there are several factors 
which throw into question whether the K-band image is a 
good constant mass-to-light ratio map of the stellar mass 
distribution despite the color correction we applied. How- 
ever, the two major factors - dust extinction and the pop- 
ulation of red super giant stars - tend to affect the arm- 
intcrarm contrast, rather than the location of the K- light 
spiral arms. Dust lanes lie preferentially inside the m = 2 
component of a galaxy's spiral (Grosb0l ct al. 2000) ab- 
sorbing interarm light, while the red super giants may ac- 
tually have their highest density in the spiral arms directly, 
where they emerged from the fastest evolving OB stars. 
This effect should become apparent in the simulations as 
slightly wrong amplitudes of the gas wiggles, leaving their 
radial position mostly unchanged. So even if the K-band 
images might actually include unaccounted mass-to-light 
ratio variations, they most probably introduce only small 
errors, which should not result in an overall mismatch of 
the models with the data. The color-corrected K-band 
images should therefore reflect the stellar mass accurately 
enough for our analysis. Certainly, there appears to be no 
better practical mass estimate that we could use for our 
analysis. 

5.2. Are there Systematic Errors in the Modelling? 

The most critical part of this study are the several mod- 
elling steps required to predict the gas velocity field for 
the comparison with the data. We apply a spatial filter 
to the K-band image before calculating the disk potential 
to reduce the significance of the clumpy H II regions (as 
described in Sect. 3.1). The residual map of the discarded 
component shows all the bright H II regions in the disk, 
giving us confidence that we have excluded much of the 
structures that reflect small-scale M/L variations. This 
correction removes roughly 3.5 % of the total K-band light 
and does not depend strongly on the number of Fourier 
components used for the fit. For two extreme decomposi- 
tions employing 6 and 16 Fourier components, the mean 
resulting relative discrepancy in the derived potentials is 



only « 10~ 6 . One remaining concern in calculating the 
potential might be that the mass density map cuts at the 
border of the image. However, the galaxy completely fits 
in the frame, fading into noise before the image border. 
Moreover, we do not perform the simulations for the com- 
plete galaxy, but only for the inner 11.6 kpc in radius. So 
the edge cutoff effect is very small and affects the part of 
the potential we arc looking at even less. 

For the dark matter component we chose an isothermal 
halo with a core because of its flexibility in fitting rotation 
curves. The functional form of the dark halo profile has 
only a second order effect on the results of the simulations 
compared to the variations due to its two basic parame- 
ters R c and Woo- We decided not to distinguish between 
different dark halo profiles for the present analysis. 

The most complex step surely is the hydrodynamical 
simulation of the gas flow in the galaxy's potential. Be- 
yond the tests of the code discussed in Section 3.2, it is 
the excellent morphological agreement between the simu- 
lated gas density profiles and the observed spiral arms that 
gives credence to the results of the code. The choice of the 
grid size was mainly motivated by the desire to achieve 
reasonable computing times, and not to exceed the seeing 
resolution. To check the grid cell size's effect on the results, 
for selected cases we also ran simulations on grids with two 
times higher resolution as well as on grids with two times 
lower resolution. These simulations showed some minor 
differences. Some gas shocks had a slightly larger ampli- 
tude in the density distribution but their locations were 
essentially unchanged. The morphology and the velocity 
field of the central bar was most susceptible to the grid cell 
size. It results in small changes in the bar position angle. 
However, we anyway exclude the central region containing 
the bar from our analysis. The final x 2— fit deviation be- 
tween simulations with different grid cells ranges at about 
6% and does not change the conclusions of this paper. Ac- 
cordingly we consider it safe to perform the simulations on 
the medium resolution grid we used. 

In a real galaxy the assumption of a global constant pat- 
tern rotation speed may not be fulfilled. In particular, we 
should expect the central bar to have a different pattern 
speed. Also, the pattern might be winding slowly, rather 
than being fixed in a corotating frame. Furthermore, the 
gas may not have a uniform temperature. If these simpli- 
fications were relaxed, the location of the spiral arms in 
the hydro-simulations would change, and eventually lead 
to a different overall fit quality. Lacking any solid basis to 
constrain these parameters, we are unable to implement 
these effects into the modelling procedure. Finally, the 
code does not include gas self-gravity The effect of gas 
self-gravity is difficult to quantify without actually per- 
forming simulations, but from the literature we know that 
it tends to amplify the gas response. Gas self-gravity also 
suppresses the tendency of the gas to shock (Lubow et al. 
1986). Since we are interested in the strength of the gas 
response to the gravitational potential and we already find 
that for high disk M/L the response is too strong, we as- 
sume that our upper limit holds also if gas self-gravity was 
included. Given the good morphological match, we may 
confidently assume that all our approximations are not far 
off. 

This leaves non-gravitationally induced gas motions as 
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the main contamination in the comparison of the mod- 
els to the observed kinematics. The influence of non- 
gravitational gas motion is in our case increased by the 
circumstance that we observe NGC 4254 from a perspec- 
tive which is not far from face-on, so any motion along the 
z-direction is expected to appear clearly in the spectra. 
Eliminating these wiggles from the rotation curve is not 
possible since we have no means to reliably identify them. 
Any method we apply for excluding parts of the data will 
be affected by some kind of bias. 

5.3. Is the Galaxy Suited for this Analysis? 

Finally we must consider the possibility that the galaxy 
we picked for our analysis might not be as suited as it 
appeared to be. NGC 4254 is not the prototype of a 
classical grand design spiral galaxy in optical wavelengths. 
However, as it can be seen in the central frame of Fig. 5, 
the galaxy exhibits in the M/L corrected K'-band image 
mainly an m = 2 spiral pattern with a strong symmetric 
part, that ends at « 5.5 kpc and fainter outer extensions. 
With its large angular size and moderate inclination, NGC 
4254 seems to be one of the most promising candidates for 
this kind of study in our sample. However, as discussed 
earlier in Section 2.2, there are indications that this galaxy 
might not be as isolated and undisturbed as one might ex- 
pect. In fact, the morphology itself implies some pcrturba- 
tive event in its recent evolution history: NGC 4254 shows 
a clear m = 1 mode and a lopsided disk. In this respect 
it was argued earlier that in-falling H I gas clumps, which 
are visible in radio data and do not emit in Ha, might be 
responsible for triggering deviations from pure grand de- 
sign structure (Phookun et al. 1993). So have we reason 
to believe that NGC 4254 is far from equilibrium? This 
is hard to tell, because on the other hand we find plenty 
of arguments that a stable propagating density wave in 
NGC 4254 is responsible for its morphology. This galaxy 
shows many similarities to the spiral galaxy NGC 5247 
whose morphological and dynamical properties were dis- 
cussed by Patsis et al. (1997), based on SPH simulations. 
We note that in our best model, the strong bisymmetric 
part of the K' spiral, ends well inside the corotation radius, 
although fainter extensions reach out to it. This picture 
is in agreement with the 4:1 SPH models of Patsis et al. 
(1997). Assuming corotation close to the characteristic bi- 
furcation of the arms at « 5 kpc on the other hand, we do 
not obtain satisfactory results (upper left frame in Fig. 5). 
Based on this is seems appropriate to conclude that NGC 
4254 is at least close to an equilibrium state and suited for 
a case study. 

6. CONCLUSIONS 

We performed hydrodynamical simulations to predict 
the gas velocity field in a variety of potentials for the spi- 
ral galaxy NGC 4254 and compared them to observations. 
These potentials consisted of different combinations of lu- 
minous (non-axisymmetric spiral) and dark matter (ax- 
isymmetric) components. The resulting gas spiral mor- 
phology reflects very accurately the morphology of the 
galaxy and allows us to specify the corotation radius or 
the pattern speed of the spiral structure quite precisely. It 
is noteworthy that within the error range given, the best 
matching pattern speed does not depend on the mass frac- 



tion of the stellar disk relative to the dark halo. For NGC 
4254 we find that the corotation lies at 7.5 ± 1.1 kpc, or 
at about 2.1 exponential K'-band disk scale lengths. From 
the kinematics of the gas simulations we could rule out 
a maximal disk solution for NGC 4254. Within the half- 
light radius the dark matter halo still has a non-negligible 
influence on the dynamics of NGC 4254: specifically, our 
fraction fd < 0.85 implies that > 1/3 of the total mass 
within 2.2 K-band disk scale lengths is dark. However, 
the comparison of the simulated gas velocity field to the 
observed rotation curves turned out to be a delicate mat- 
ter. The observed rotation curves show a significant num- 
ber of bumps and wiggles, presumably resulting from non- 
gravitational gas effects, that complicate the identification 
of wiggles induced by the massive spiral arms. Therefore, 
beyond concluding that the disk is less than 85% of max- 
imal, we were unable to specify a particular value for the 
disk mass or to test the results from Bottema (1997) or 
Courteau and Rix (1999). But already with this statement 
we differ from the conclusions of Debattista and Scllwood 
(2000) and Weiner et al. (2001b), who argue that their 
conclusions for maximal disks of barred galaxies also hold 
for non-barred spirals. Since we only analyzed one galaxy 
so far it is inappropriate to state here that the centers of 
unbarred spirals, what we still consider NGC 4254 to be 
despite its small bar-like structure in the very center, are 
generally governed by dark matter. 

In the near future we will extend this analysis to more 
galaxies from our sample. This will put us into the position 
to decide, whether single galaxies differ very much in their 
dark matter content, or if the bulk of the spirals show 
similar characteristics. 

We thank Panos Patsis and Julien Devriendt for their 
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